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Abstract: We present an analysis of the global neutrino oscillation data including the 
LSND result in terms of four-neutrino mass schemes and a CPT violating three-neutrino 

r- -h 

' framework. We find that the strong preference of oscillations into active neutrinos im- 

plied by recent solar+KamLAND as well as atmospheric neutrino data allows to rule out 
£SJ ■ (2+2) mass schemes at 5.8er, whereas (3+1) schemes are disfavoured at the 3.2a level by 

short-baseline experiments. The CPT violating scenario is disfavoured at 3.4er by Kam- 
LAND and atmospheric anti-neutrino data. We conclude that currently no convincing 
explanation for the LSND signal exists, and it remains a puzzle how to reconcile it with 
*^ the other evidences for neutrino oscillations. 



1. Introduction 

The large amount of data from the solar |l|, ^] , atmospheric |3], Q , KamLAND reactor || and 
K2K accelerator |6] neutrino experiments is beautifully described in terms of three-flavour 
neutrino oscillations 0]. To reconcile also the evidence for oscillations from the LSND ex- 
periment [|| with all the other data provides a big challenge for neutrino phenomenology. 
To obtain the required three neutrino mass-squared differences of different orders of magni- 
tude it has been proposed to introduce a light sterile neutrino H, or to allow for violation 



of CPT invariance by neutrino masses and mixing |1C], 11 1. Here we present an analysis of 
the global neutrino oscillation data in terms of four-neutrino and CPT violating schemes, 
including data from solar, atmospheric, KamLAND and K2K neutrino experiments, the 
LSND experiment, as well as data from short-baseline (SBL) experiments 13, 14 1 and 
reactor experiments |TH] reporting no evidence for oscillations. 



After fixing some notations in Sec. |2|we show in Sees. ||, ||, |5]that for all possible types of 
four-neutrino schemes different sub-sets of the data are in serious disagreement and hence, 
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Figure 1: The six four-neutrino mass spectra, divided into the classes (3+1) and (2+2). 



four-neutrino oscillations do not provide a satisfactory description of the global oscillation 



data including LSND. The details of our calculations can be found in Refs. 16, 17, 18]. In 
Sec. we show that also the CPT violating scenario is strongly disfavoured [19|. In Sec. ^ 
we briefly comment on some proposed alternative solutions for the LSND puzzle, and we 
summarize our results in Sec. || 



2. Four-neutrino schemes: notations and approximations 

Four-neutrino mass schemes are usually divided into the two classes (3+1) and (2+2), as 
illustrated in Fig. [I]. We note that (3+1) mass spectra include the three-active neutrino 
scenario as limiting case. In this case solar and atmospheric neutrino oscillations are 
explained by active neutrino oscillations, with mass-squared differences Am^ 0L and Am^ TM , 
and the fourth neutrino state gets completely decoupled. We will refer to such limiting 
scenario as (3+0). In contrast, the (2+2) spectrum is intrinsically different, as there must 
be a significant contribution of the sterile neutrino either in solar or in atmospheric neutrino 
oscillations or in both. 

Neglecting CP violation, in general neutrino oscillations in four-neutrino schemes are 
described by 9 parameters: 3 mass-squared differences and 6 mixing angles in the unitary 
lepton mixing matrix. Here we use a parameterisation introduced in Ref. [p|, which 
is based on physically relevant quantities: the 6 parameters Amg 0L , 6 SOL , Am^ TM , 6 ATM , 
^ m LSND' #lsnd are similar to the two-neutrino mass-squared differences and mixing angles 
and are directly related to the oscillations in solar, atmospheric and the LSND experiments. 
For the remaining 3 parameters we use ri s ,rj e and d^. Here, n s (r/ e ) is the fraction of 
u s (y e ) participating in solar oscillations, and (1 — d^) is the fraction of participating in 
oscillations with Am^ TM (for exact definitions see Ref. [fL6| ). For the analysis we adopt the 
following approximations: 



1. We make use of the hierarchy Amg 0L <C Amf XM <C Am^ SND . This means that for 
each data set we consider only one mass-squared difference, the other two are set 
either to zero or to infinity. 
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Figure 2: Parameter dependence of the different data sets in our parameterisation. 



2. In the analyses of solar and atmospheric data (but not for SBL data) we set r] e = 1, 
which is justified because of strong constraints from reactor experiments [14, [la ]. 



Due to these approximations the parameter structure of the four-neutrino analysis 
gets rather simple. The parameter dependence of the four data sets solar, atmospheric, 
LSND and NEV is illustrated in Fig. ^. The NEV data set contains the experiments KAR- 
MEN CDHS Hi, Bugey @ and CHOOZ/Palo Verde p|, reporting no evidence for 
oscillations. We see that only rj s links solar and atmospheric data and d^ links atmospheric 
and NEV data. LSND and NEV data are coupled by Attend and $lsnd- 



3. (2+2): ruled out by solar and atmospheric data 

The strong preference of oscillations into active neutrinos in solar and atmospheric oscilla- 
tions leads to a direct conflict in (2+2) oscillation schemes. We will now show that thanks 
to recent solar neutrino data (in particular from SNO ||) in combination with the Kam- 
LAND experiment 0, and the latest SK data on atmospheric neutrinos Q the tension in 
the data has become so strong that (2+2) oscillation schemes are essentially ruled out. 1 

In the left panel of Fig. || we show the A% 2 from solar neutrino data as a function of rj s , 
the parameter describing the fraction of the sterile neutrino participating in solar neutrino 
oscillations. One can see that the improved determination of the neutral current event rate 
from the solar 8 B flux implied by the salt enhanced measurement in SNO leads to a 
significant tightening of the constraint on r/ s : the 99% CL bound improves from rj s < 0.44 
for pre-SNO salt to rj s < 0.31 including also the most recent SNO data. Moreover, large 
values of rj s > 0.5 are now ruled out by many standard deviations. The outstanding results 
of the KamLAND reactor experiment Q confirmed the LMA solution of the solar neutrino 



problem [2C, 22]. Apart from this very important result the impact of KamLAND for an 
admixture of a sterile neutrino is rather limited (see Fig. |3|). Since KamLAND on its own 
has no sensitivity to rj s the bound is only indirectly affected due to the better determination 
of Amg 0L and 9 S0L . The combined analysis leads to the 99% CL bound rj s < 0.27. 



1 Details of our analyses of the solar, KamLAND and atmospheric neutrino data can be found in Refs. |lg| , 
fed]. For an earlier four-neutrino analysis of solar and atmospheric data see Ref. ]21|] . 
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In contrast, in (2+2) schemes atmospheric data imply r/ s > 0.65 at 99% CL, in clear 
disagreement with the bound from solar data. In the right panel of Fig. |3] we show the 
Ax 2 for solar data and for atmospheric combined with SBL data as a function of rj s . 
Furthermore, we show the x 2 °f the global data defined by 

X 2 (rjs) = A X 2 S0L (Vs) + Axi TM+S BL(%) • (3-1) 

From the figure we find that only at the 4.7a level a value of rj s exists, which is contained 
in the allowed regions of both data sets. This follows from the x 2_vaiUe Xpc = ^2.4 



shown in the figure. In Refs. [17, |23| we have proposed a statistical method to evaluate 
the disagreement of different data sets in global analyses. The parameter goodness of fit 
(PG) makes use of the x 2 defined in Eq. ( |3.1[ ). This criterion evaluates the GOF of the 
combination of data sets, without being diluted by a large number of data points, as it 
happens for the usual GOF criterion (for details see Ref. f23|). We find Xpc = Xmin = 33.2, 
which corresponds to 5.7a. We conclude that (2+2) mass schemes are ruled out by the 
disagreement between the latest solar and atmospheric neutrino data. This is a very robust 
result, independent of whether LSND is confirmed or disproved. 2 



4. (3+1): strongly disfavoured by SBL data 

It is known for a long time [25] that (3+1) mass schemes are disfavoured by the comparison 



of SBL disappearance data pL |l|] with the LSND result. In Ref. p§ we have calculated 



2 Sub-leading effects beyond the approximations adopted here should not affect this result significantly. 
Allowing for additional parameters to vary at the percent level might change the ratio of some observ- 
ables (53], however, we expect that the absolute number of events relevant for the fit will not change 
substantially. 
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Figure 4: Upper bound on sin 2 29 LSND from NEV and atmospheric neutrino data in (3+1) 
schemes |2(| compared to the allowed region from global LSND data || and decay-at-rest (DAR) 
LSND data §7|. 



an upper bound on the LSND oscillation amplitude sin 2 2# LSND resulting from NEV and 
atmospheric neutrino data. From Fig. |I] we see that this bound is incompatible with the 
signal observed in LSND at the 95% CL. Only marginal overlap regions exist between 
the bound and global LSND data if both are taken at 99% CL. An analysis in terms 
of the parameter goodness of fit [fl?]] shows that for most values of Am^ SND NEV and 
atmospheric data are compatible with LSND only at more than 3a, with one exception 
around Am^ SND ~ 6 eV 2 , where the PG reaches 1%. These results show that (3+1) schemes 
are strongly disfavoured by SBL disappearance data. 



5. Comparing (3+1), (2+2) and (3+0) hypotheses 



With the methods developed in Ref. |16[ we are able to perform a global fit to the oscillation 
data in the four-neutrino framework. This approach allows to statistically compare the 
different hypotheses. Let us first evaluate the GOF of (3+1) and (2+2) spectra with the 



help of the PG method described in Ref. [23|. We divide the global oscillation data in the 



4 data sets SOL, ATM, LSND and NEV. Following Ref. g7| we consider 

X 2 = Ax 2 ol(#sol, Am 2 0L , T) s ) + Ax! tm (0atm, Am 2 TM , <r) s , dp) 
+ AxLv(<W>, Am 2 SND , rj e ) + Ax 2 SND (6> LSND , Am 2 SND ) , 



(5.1) 
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SOL ATM LSND NEV 


Xpg 


PG 


(3+1) 
(2+2) 


0.0 0.4 7.2 7.0 
4.1 29.2 1.2 9.7 


14.6 
44.1 


5.6 x 10" 3 (3.2a) 
6.1 x 10" 9 (5.8a) 



Table 1: Parameter GOF and the contributions of different data sets to Xpg m (3+1) an d (2+2) 
neutrino mass schemes. 



where A X 2 X = Xx - (Xx)min (X = SOL, ATM, NEV, LSND). In Tab. g we show the 
contributions of the 4 data sets to Xpg = Xmin f° r (3+1) and (2+2) oscillation schemes. As 
expected we observe that in (3+1) schemes the main contribution comes from SBL data 
due to the tension between LSND and NEV data in these schemes. For (2+2) oscillation 
schemes a large part of Xpg comes from solar and atmospheric data, however, also SBL 
data contributes significantly. This comes mainly from the tension between LSND and 
KARMEN |p7fl , which does not depend on the mass scheme and, hence, also contributes 
in the case of (2+2). Therefore, the values of Xpg m Tab. || for (2+2) schemes are higher 
than the one given in Sec. pi where the tension in SBL data is not included. 

The parameter goodness of fit is now obtained by evaluating Xpg f° r ^ DOF [23]. 
This number of degrees of freedom corresponds to the 4 parameters n s ,d^,9 LSNB , Am 2 SND 
describing the coupling of the different data sets (see Eq. (|5.1|) and Fig. ||). The best GOF 
is obtained in the (3+1) case. However, even in this best case the PG is only 0.56%. The 
PG of 6.1 x 10 -9 for (2+2) schemes shows that these mass schemes are essentially ruled 
out by the disagreement between the individual data sets. 

Although we have seen that none of the four-neutrino mass schemes can provide a 
reasonable good fit to the global oscillation data including LSND, it might be interesting 
to consider the relative status of the three hypotheses (3+1), (2+2) and the three-active 
neutrino scenario (3+0). This can be done by comparing the x 2 values of the best fit point 
- which occurs for the (3+1) scheme - with the one corresponding to (2+2) and (3+0). 
First we observe that (2+2) schemes are strongly disfavoured with respect to (3+1) with 
a Ax 2 = 29.5. For 4 DOF this is equivalent to an exculsion at 4.5o\ Furthermore, we 
find that (3+0) is disfavoured with a Ax 2 = 20.0 (99.95% CL for 4 DOF) with respect to 
(3+1). This reflects the high statistical significance of the LSND result, since in a (3+0) 
no effect is predicted for LSND. 



6. CPT violation 



The CPT violating solution for the LSND results proposed in Refs. [ 10 1 is based on the 
observation that the signal in favour of oscillations in LSND is dominated by the anti- 
neutrino signal from decay-at-rest data. The neutrino data from the decay-in-flight event 
sample gives a — > v e transition probability of (0.10 ± 0.16 ± 0.04)% which, although 
consistent with the anti-neutrino signal of Pp M ^p e = (0.264 ± 0.067 ± 0.045)%, is also 
perfectly consistent with the absence of neutrino oscillations in LSND. Therefore, it has 
been proposed to adopt different masses and mixings for neutrinos and anti-neutrinos 
to reconcile the LSND result (n], 28] without introducing a sterile neutrino. In such a 
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Figure 5: CPT violating neutrino mass scheme proposed in Ref. 

scheme solar and atmospheric neutrino oscillations are driven by Am 2 0L and Am^ TM , and 
atmospheric and LSND anti-neutrino oscillations by Am^ TM and Ara^ SND . 

This scenario is strongly affected by the KamLAND results on anti-neutrino disappear- 
ance, which in combination with the constraints from CHOOZ indicate oscillations with 
Am 2 < 10~ 4 eV 2 . Since this requires a third mass scale in the anti-neutrino sector the CPT 



violating solution seems to be disfavoured. However, in Ref. 11] it was claimed that the 
mass spectrum shown in Fig. | still may provide a fit to the data, although no anti-neutrino 
mass-squared difference in the range required for atmospheric oscillations is present. This 
claim is based on the fact that Super-Kamiokande cannot distinguish between neutrinos 
and anti- neutrinos, and the atmospheric data sample is dominated by neutrinos by a factor 
2-4 over anti-neutrinos. 



In Ref. 1 19 ] we have performed an anlysis of the global oscillation data in the framework 
of the CPT violating neutrino scheme. First we have analized all data except LSND - 
from solar, atmospheric, reactor (KamLAND, CHOOZ, Bugey) and K2K experiments - 
in a general three-flavour context, by allowing for different neutrino and anti- neutrino 
oscillation parameters. We found that there is no evidence for any CPT violation in that 
data set. In fact, the best fit point occurs very close to CPT conservation, and the best fit 
in the CPT conserving case has Ax 2 of only 0.5. 

In Fig. | we show the allowed regions of the all-but-LSND data compared to the LSND 
allowed regions, projected onto the plane of the LSND oscillation parameters Am 2 SND = 
Amjj and sin 2 2# LS nd = sin 2 #23 sin 2 26*13. This figure illustrates that below 3a CL there 
is no overlap between the allowed region of the LSND analysis and the all-but-LSND one, 
and that for this last one the region is restricted to Am 2 SND < 0.02 eV 2 . At higher CL 
values of Am 2 gND ~ C(eV 2 ) become allowed - as determined mainly by the constraints 
from Bugey - and an agreement becomes possible. We find that in the neighbourhood 
of Am 2 SND = 0.9 eV 2 and sin 2 2(9 LSND = 0.01 the LSND and the all-but-LSND allowed 
regions start having some marginal agreement slightly above 3a CL (at Ax 2 = 12.2). A 



-7- 



International Workshop on Astroparticle and High Energy Physics 



Thomas Schwetz 




sir,22e LSND = sin 2 e 23 sin 2 2e i3 



Figure 6: 90%, 95%, 99%, and 3cr CL allowed regions (filled) in the (At 



sin 2 20, 



plane 



required to explain the LSND signal together with the corresponding allowed regions from our 
global analysis of all-but-LSND data. The contour lines correspond to A\ 2 — 13 and 16 (3.2<r and 
3.6cr, respectively). 



less fine-tuned agreement appears at 3.3cr CL (A X 2 ~ 14) for Am LSND > 0.5 eV 2 and 
sin 2 26 LSND < 0.01. 

Alternatively the quality of the joint description of LSND and all the other data can be 



evaluated by performing a global fit based on the total x 2 -f unc ti° n 



Xall-but-LSND + 



XlsnD' an d applying a goodness-of-fit test. The best fit point of the global analysis is 
sin 2 2# LSND = 6.3 x 10~ 3 and Am 2 SND = 0.89 eV 2 . In the following we will again use the 
parameter goodness-of-fit [ )23[ |, which is particularly suitable to test the compatibility of 
independent data sets. Applying this method to the present case we consider the statistic 



X 



—but— LSND 

(b.f.) + A x 2 SND (b.f.), 



(6.1) 



where b.f. denotes the global best fit point. The x 2 °f Eq. ( |6.1| ) has to be evaluated for 2 
dof, corresponding to the two parameters sin 2 2# LSND and Arajj = Am 2 SND coupling the two 
data sets. From Ax 2 u_ but _ L gND = 12.7 and Ax 2 SND = 1.7 we obtain x 2 = 14-4 leading to 
the marginal parameter goodness-of-fit of 7.5 x 10~ 4 , corresponding to 3.4<r. 

These results show that atmospheric data is sensitive enough to anti-neutrino oscilla- 
tions, and hence, no reasonable fit can be obtained without a mass-squared difference of the 
order 10 -3 eV 2 for anti- neutrinos. Consequently the CPT violation scenario is disfavoured 
by global anti-neutrino data from reactor, atmospheric and the LSND experiments. 
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7. Alternative solutions 

In this section we briefly mention some alternative solutions which have been proposed to 
explain the result from LSND. 



In Ref. [29 1 a lepton number violating decay of the muon has been considered as the 
source for the LSND signal. In this model the v e event excess in LSND stems from 
the decay fi + — > e + + v e + D x , where v x can be an anti- neutrino of any flavour. This 
scenario seems to be in disagreement with the non-observation of any v e excess in 
KARMEN. 

Recently an explanation in terms of five neutrinos has been proposed |i(J]. In such a 
(3+2) scheme solar and atmospheric data are explained dominantly by active oscil- 
lations, similar to the (3+1) case. The fit of SBL data is improved with respect to 
(3+1) by involving two large mass splittings Am^ ~ 0.9 eV 2 and Am^ ~ 20 eV 2 . 
In this case it seems to be very dificult to reconcile these large neutrino masses with 



constraints from cosmology [31] 



In Ref. J32J it has been pointed out that involving a sterile neutrino and CPT violation 
simultaneously all the data can be explained. Once such drastical modifications of 
standard physics are accepted, one may adopt a (3+1) like mass scheme, even with 
identical masses for neutrinos and anti-neutrinos, but allowing for different mixings. 
Since CDHS constrains only SBL muon neutrino disappearance (not anti- neutrinos) , 
whereas Bugey restricts v e mixing, a small violation of CPT would suffice to get 
around the SBL bounds and explain the LSND signal. 



8. Summary 

Performing a global analysis of current neutrino oscillation data we find that neither for 
four-neutrino schemes nor for CPT violation in a three-neutrino framework a satisfactory 
fit to the data is obtained. 

• The strong rejection of non-active oscillation in the solar +KamL AND and atmo- 
spheric neutrino data rules out (2+2) schemes, independent of whether LSND is 
confirmed or not. Using an improved goodness of fit method especially sensitive to 
the combination of data sets we find that (2+2) schemes are ruled out at the 5. So- 
level. 

• (3+1) spectra are disfavoured by the disagreement of LSND with short-baseline dis- 
appearance data, leading to a marginal GOF of 5.6 x 10~ 3 (3.2a). If LSND should 
be confirmed it would be very desirable to have more data on u e and/or SBL 
disappearance to decide about the status of (3+1). In that case a positive signal is 
predicted right at the sensitivity edge of existing experiments. 

• Also in the case of the CPT violating three-neutrino mass scheme no satisfactory fit 
can be obtained. This scheme is disfavoured by KamLAND, CHOOZ, Bugey and 
atmospheric anti-neutrino data at the 3. 4a level. 
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Alternative solutions, which have been proposed to reconcile the LSND result, seem 
to be in disagreement with some data, or require a very radical modification of standard 
physics. We conclude that currently no convincing explanation for the LSND result exists, 
and it remains a puzzle how to reconcile this evidence with the rest of the global data. 
Therefore, it is very important to settle this issue experimentally. A confirmation of the 
LSND signal by the MiniBooNE experiment [33] would be very exciting and would require 
some new ideas. 
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